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Abstract. An analytic model a la Middleton of the impulsive noise component in the data 
of interferometric gravitational wave detectors is proposed, based on an atomic representation 
of glitches. A fully analytic characterization of the coherent network data analysis algorithm 
proposed by Rakhmanov and Klimenko is obtained, for the simplest relevant case of triggered 
detection of unmodeled gravitational wave bursts, using the above noise model. The detector's 
performance is evaluated under a suitable central-limit hypothesis, and the effects of both the 
noisiness of the pseudo-templates, and the presence of the impulsive noise component are 
highlighted. 



1. Introduction 

Gravitational wave (henceforth GW) astronomy is expected to open an essentially new 
observational window on the physical Universe. Several classes of GWs of cosmic origin are 
currently being sought for, including continuous, transient and stochastic ones. An essential 
distinction among these different signals concerns our ability in modeling the expected 
waveforms. GW bursts (henceforth GWB) are a paradigm of transient signals for which 
only a few physically -based models exist Q-G). 

GW detectors (with specific reference to present-day large baseline optical interferom- 
eters) are invariably affected by transient disturbances of various origin [4 ]. Using auxiliary 
channels to monitor the status of the instrument and its environment may help identifying and 
vetoing these disturbances. Experimental evidence suggests that a residual impulsive com- 
ponent will nonetheless be present in the data. Distinguishing these spurious noise glitches 
from true GWB of cosmic origin will be almost impossible, when only data from a single 
detector are available. It becomes feasible, in principle, if the outputs of several detectors are 
suitably combined. Using data from several detectors it is further possible to reconstruct the 
GW signal waveform, which encodes the relevant source physics, so as to capitalize on and 
refine astrophysical models. 

As an historical heritage of acoustic GW detectors 0, various coincidence algorithms, 
based on consistency tests among candidate-events gathered by different detectors, have been 
studied and tested [|6|-|l7j. These algorithms, while conceptually simple and computationally 
inexpensive, turn out to be less efficient, in general, compared to coherent techniques, 
where the output data from several sensors are combined to form a suitable detection 
statistic to be used in classical hypotheses tests [8j. Several coherent techniques have been 
hitherto proposed Il9l-lfl6l. but only a few (e.g., WAVEBURST flj], X-PIPELINE fl3] and 
RIDGE lfl6ll ) have been fully implemented to date in the data analysis pipelines of running 
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experiments. Considerable work is still needed for completely characterizing alternative 
coherent algorithms in terms of performance and computational cost. 

This paper (the first in a suite, where we propose to investigate problems of increasing 
complexity) attempts to provide a quantitative answer to the rather fundamental question 
of how well a network of several GW detectors may discriminate true GWBs from local 
disturbances (glitches) using coherent detection statistics. Among the essential benefits 
provided by coherent network operation, we mention: i) the ability of detecting unmodeled 
signals; ii) the capability of rejecting local disturbances; the possibility of iii) retrieving the 
source position on the celestial spere and iv) reconstructing the gravitational waveforms. Here 
we focus on the first two properties, assuming for simplicity that the GW direction of arrival 
(henceforth DOA), and the time of occurrence (henceforth TOO) of the event are known 
(triggered search) from observations of different nature (e.g., electromagnetic, neutrino, etc.). 

To this end, modeling the impulsive noise component (the glitches) is a key, and yet 
open, issue. In this paper we adopt, for the first time to the best of our knowledge, a general 
representation of impulsive noise proposed by D. Middleton in a series of seminal papers 
ifTTll . Glitches are accordingly modeled as time-frequency atoms, i.e., transients whose energy 
content is almost confined to a compact region in the time frequency plane, and characterized 
in terms of a few (random) parameters; we adopt the possibly simplest (though observation- 
driven) model for such atoms: real-valued sine-gaussian (SG) functions. The impulsive noise 
component in each detector, is modeled as a random train of these atoms, each occurring 
independently on the other ones. 

In order to keep the formal complexity to a minimum, while still capturing the key 
aspects of the problem, we make a number of simplifying assumptions, summarized below. 

We refer to the coherent cross-correlation method introduced by Rakhmanov and 
Klimenko (henceforth RK) in [ 18 1 and deduce its statistical properties in analytic form. This 
latter may be recognized as a natural extension of the matched correlator concept to the 
detection of unmodeled GWBs with a redundant network of detectors. 

We limit to the simplest (though realistic and up-to-date) case of a network composed of 
three interferometers of comparable sensitivity, with specific reference to the LIGO-Hanford 
(henceforth LH), LIGO-Livingston (henceforth LL) and Virgo (henceforth V) detectors, and 
restrict to the case where the incident wave is linearly polarized. 

This paper is accordingly laid out as follows. In Section 2 we recall the RK formalism, 
and deduce the distribution of its detection statistics under the H\ hypothesis (GWB in the 
data). In Section 3 the adopted atomic representation of glitches is briefly introduced, and the 
first two moments of the RK detection statistics under the Hq hypothesis (only noise in the 
data) are derived following a simple heuristic reasoning, making the simplifying assumption 
that in each detector no more than a single glitch may occur in the analysis window. This 
assumption is relaxed in Section 4, where we propose the anticipated fairly general and 
rigorous approach a la Middleton to model the impulsive component of the interferometer 
noise; the results obtained in Sections 3 and 4 are shown to coincide, under the appropriate 
simplifying assumptions. In Section 5, based on the above results, we evaluate the RK 
correlator based detector's performance by numerical experiments. Conclusions and hints 
for future work follow under Section 6. 

2. The RK Coherent Analysis Algorithm 

Whenever the sought waveform is known a priori, optimal detection in additive stationary 
(band-limited) white gaussian noise is achieved by matched-filtering the data with a template 
of the sought waveform. The output of the filter (also known as matched correlator) has to 
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be compared to a properly chosen threshold in order to decide about the presence or absence 
of the signal in the data. When the signal shape is known, except for a finite number of 
parameters, a set of correlators corresponding to a suitably dense covering of the parameter 
space can be computed, and the largest one exceeding the threshold selected, yielding an 
estimate for the signal parameters. 

For unmodeled GWBs the matched filtering technique cannot be adopted. It is thus 
basically impossible to distinguish a GWB from a spurious glitch surviving the auxiliary- 
channel based vetos, in the data of a single interferometer. 

One possible way to circumvent this difficulty using data from a (redundanfl) network 
of detectors has been proposed in [18], and will be shortly recalled hereinafter. 

Let a (plane) gravitational wave with linearly polarized (TT gauge) components h + (t) 
and h x (t) impinge on a set of interferometric detectors located at r = r,-, i = 1,2,3, from a 
direction £l s . 

In the absence of noise, the detector outputs can be writterU 

Sj(0 = F?(a s )h+[t - Ufa)] +F? (0.)A X [* ~ *<(^)L i = 1-2,3, (1) 

where T ( (n s ) = c~ n ■ r; is the propagation delay of the (plane) wavefront referred to its arrival 
at the spatial origin (usually taken coincident with the Earth center), n is the unit wave-vector, 
and Fi~(Q. s ), F* (£l s ) are the pattern functions describing the antenna directional response. 
Equations (Q~|i can be rewritten in matrix form as follows, 

, / u . (t\ \ 

(2) 

V S 3 (t) J V F 3 + (^) F 3^) J V " Xl ' ; y 
where 

S,(t) = S' i [t + x i (a s )} 1 i= 1,2,3 (3) 

are the properly time-shifted noise-free detector outputs, and the matrix on the r.h.s. is the 
network response matrix. The rank of this latter cannot exceed 2, hence, 

St F+ F* 

det | S 2 F 2 + F 2 x I = 0. (4) 

S3 F+ F x 

Expanding the determinant in the elements of the first column, we obtain the null condition 

m 




where 



A l (£l s )S l (t)+A 2 (£ls)S2{t)+A 3 (£l s )S 3 (t) =0, (5) 
At =F+F 3 X -F+F* 

A 2 =F+F* -F+F*. (6) 
A 3 = FfFf - F 2 + F l x 

In (O (and whenever possible, hereafter) the dependence on Q. s of A, and F^' x is omitted for 
notational ease. 

$ A network of detectors is redundant if, in the absence of noise, the output of each detector can be expressed in 
terms of the outputs of the others. Since a gravitational signal has only two independent polarization components, 
any network of three or more (differently located and oriented) detectors is redundant in the common observational 
band. 

§ In writing eq. {TJ, we make the usual (competing) assumptions that the gravitational wave signal is i) short enough 
to ignore the reciprocal motion between source and detector negligible, and ii) spectrally narrow enough to make the 
response of the antenna practically instantaneous. 
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From eq. Q one may infer that, in the absence of noise, the output of each and any detector 
in the network is proportional to a linear combination of the outputs of the remaining two, 
namely 

Li(t) = -AjSj(t) -AtS t (t) =A i S i (t), i,j,£= 1,2,3, j?£? i. (7) 

and is thus a templat^ for S 1 ,-. The actual interferometer outputs V,(? ), however, differ from 
the Si(t) due to the presence of noise, viz. 

V i (t)=S l (t)+n i {t), i = 1,2,3. (8) 

Accordingly, by using the Vj in place of 5, in (0 we obtain a noisy template for S 1 ,-, 

Wi{t) = -AjVj{t)-A e V e (t) =Ai [Si(t) + Vi(t)} , 

ij,£=l,2,3,j?£?i, (9) 

where 

v i (t) = ~\Ajn J (t)+A l ni(t)], i,j,£= 1,2,3, i. (10) 

The Wj(t) can thus be used, following RK lfl8l . to compute the pseudo matched-correlators 

Ci=(Vi,Wi) = [ Vi(t)Wi(t)dt « 
J©, 

£ W^/r'VrWf, i= 1,2,3. (11) 

m= 1 

In (fTTT> ©, is the (T seconds wide) analysis window, f s the sampling frequency, A( s = [f s T\ 
the number of samples in ©,-, V, m and W lm the time samples of V (t) and W(f), respectively, 
V, = {Vn,V i2 , . . .,ViN s }, and W,- = {W n ,W i2 , . . . ,W iNs }. 

In view of the relatively large value of N s in (fTTT i. typically ~ 10 2 , in the following 
we shall make the working assumption that a suitable form of the (generalized) Central 
Limit Theorem may be invoked ll20l to argue that the C, are normal distributed under both 
hypotheses H\ (signal present) and Hq (no signal), despite the presence of glitches, which 
makes the interferometer noises depart from Gaussianity. Only the first two moments will be 
accordingly needed to characterize them. 

2.1. RK Correlators Distributions under Hy 

In this Section we assume that the occurrente of a GWB and an instrumental glitch in the 
same analysis window can be neglected, being extremely unlikeljU Denoting the moments 
for the Hi case with the subscript 1, one readily obtains (z = 1,2,3) 

Hf ] = E [Qm = A f (0,) / Sf(t )dt « f^Ai^Si ■ Sf , (12) 
J&i 

For all practical purposes, one may use as a template for S, any quantity differing from 5, by an arbitrary 
multiplicative factor 1 19 |. 

f This (reasonable) assumption may be relaxed using the moments of the total noise derived in Sect. 4, in lieu of 
those of the Gaussian component alone. Results pertaining to this more general case will be presented elsewhere. 
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al' ] ) 2 =Var[C i \H 1 ]=Af{n s ) 



N+N 



2^2, 



Sf(t)dt- 



f s ~ 2 Aj(a s ) [(oj 2 + of) St ■ Sf + afafN s ] 



(13) 

where N and Nj denote the one-sided power spectral densities of «,•(?) and v, (f ) respectivel)I+l 
In deriving (TT2T) and (13[ we capitalize on the statistical independence between n; and v,, the 
obvious identities 



E[n i ]=E[v i ]=E[n i V i }=0, 



(14) 



and the relationship Nj = 2af/f s , valid for band-limited gaussian white noise with standard 
deviation O;. 

The performance of the RK correlator is described in terms of its deflection [19|, aka 
signal to noise ratio (SNR) defined as 



<' [(o? + a?)S r Sf + o?6?tf,.] 



1/2 : 



(15) 



In the following we shall assume for simplicity that all detectors in the network have 
comparable noise PSDs, thus letting Ni = N and a, = a, Vi. Accordingly, 



a} 



A2 



so that eq. (fT~5T > becomes 
d« = - 



A,(S,-S 



n l/2 



1/2- 



(16) 



(17) 



The deflection ( TTTb can be written in a more transparent form by introducing the 
quantities 



gW - 
°s 



N/2 



1/2 



(18) 



representing the signal to noise ratio of a perfect matched filter applied to the actual data at 
the output of detector-/, and 



U/2 



1/2 



(19) 



+ The last term in 1 1 3t is obtained from 



using the band-limited white-noise formula 

E 



I dt I dsni(t)ni(s)Vi(t)Vi(s) 

J[T] J[T] 
1 

»W»W] = I J_ g eM'2nf(t- S )]df. 
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where 

h 2 rss = J [h+{t) 2 + h x (t) 2 ]dt (20) 

is a frequently used measure of the GWB strength. For the simplest case of linearly polarized 
GWBs, 

SP = \Fi\Sk, (21) 

where f, = F, ' x , depending on the wave polarization, and 5/, represents the intrinsic signal 
to noise ratio of a perfect matched filter applied to the bare gravitational waveform embedded 
in the detector noise; such a deflection would be attained if the antenna were isotropic 
(i^ + ' x = 1), and the template noise-free. 

The deflection d^ in < fT~5T > can thus be conveniently written: 

dU=8®E i (£l S) N s ,F i ,8 h ), (22) 

where 

At 



Si(£l s ,N s ,F h 8 h ) = 



A 2 +A)+A}- 



-NsfAj+AfjQFilSh)- 



1/2- 



(23) 



measures the SNR degradation of the RK correlator C, w.r.t. the perfect matched filter acting 
on the output data of detector-i, due to the noisiness of the template. It is seen from ( 1231 that 
£;(•) depends on the DOA, the number of samples N s in the analysis window, the polarization- 
dependent pattern function Fj, and the intrinsic deflection 5/,. 

Figures la and lb display the sky maps of the function in d23b for the LH, LL and V 
detectors for the two linear polarizations, for N s = 100 and two extremal values of 5/,, namely 
5/, = 10 and 5/, = 100, respectively. The source position in Figs, la and lb is parameterized in 
terms of the polar and azimuthal angles & s , q> s in an Earth-centered coordinate system whose 
polar axis points to the North-Pole, and (p s — identifies the Prime Meridian. 



3. Glitches 



Available experimental evidence |21 1,|22| suggests that instrumental noise glitches can 
be efficiently modeled as atoms 11231 in the time-frequency plan^3 [ 24 1 . Atoms are waveforms 
with almost-compact time-frequency support. They can be characterized in terms of their 
energy content, and their first and second order moments, i.e., occurrence time to, center 
frequency /o, effective duration a t and bandwidth ay. 

The choice of an atom family (technically called a dictionary) appropriate to modeling 
glitches in GW interferometers must be compliant to and derived from experimental evidence. 
In this connection, the work in [21 1,[22|, aimed at classifying glitches and identifying glitch 
families (clusters in parameter space) is particularly relevant. It should be noted that atoms 
in general form overcomplete systems, and this fact must be taken properly into account in 
deducing the distributions of the atom parameters from observed glitch populations [25 1. 

We adopt here the possibly simplest atom, the real valued sine-Gaussian (SG) functions 
defined by 

yr(t-to;go,Jb,Oi) = gosin[27r/o(r-f )]e- ( ^ o)2 / a ' 2 , (24) 

* GWBs can also be modeled as atoms. For the detection technique adopted here, however, only the GWB energy 
is relevant, not its shape. 



Modeling the Impulsive Noise Component.. 



7 



whose waveform and time-frequency representation are shown in Figure 2. The SG atom is 
entirely characterized by its shape parameters go, /q and a t , and effective occurrence (firing) 
time fo- 

The choice of the SG dictionary is suggested by the fact that a wide variety of observed 
glitches in the data channel are well modeled as SG atoms ll22l . and is further motivated 
by its structural simplicity, minimum time-frequency spread, ((7,0/ = (4-n)~ l ), and positive - 
definiteness of its Wigner-Ville transform. These properties should likely permit to represent 
the instrumental transients in a close-to-optimal (i.e., minimally redundant) way (see, e.g., 
(261,1123). 



3.1. RK Correlators Distributions under Hq - Heuristic Approach 

In this Section we obtain a heuristic characterization of the RK correlator distribution 
under Hq by considering the glitches as (spurious) signals with random parameters. 
Specifically, we derive the "average" among the marginal distributions of the RK correlator 
corresponding to all possible glitch realizations in the network. To keep the analysis as simple 
as possible, we assume that no more than a single glitch may occur in the analysis window 
in each interferometer. This restriction will be removed in the next Section, where a fairly 
general and rigorous model for the impulsive noise component will be proposed. 

It is expedient to write the moments under Hq (denoted with a subscript 0) as follows: 

$ = (i -n) 3 jug,+n(i -nfrifi +n 2 (i -n)^ +nVg, (25) 



(a«) 2 ^(i-n)3(a«) 2 + n(i-n) 2 (a«) 2 - 



n 2 (i-n)«;r+rf a{»)", (26) 



2 



where f^\, <Jq \ refer to the cases where k detectors (k = 0,1,2,3) in the network exihibit 
a glitch within the analysis window, and the corresponding factors in front of them are the 
related occurrence probabilities, II being the (known) probability of observing a single glitch 
in the analysis window. 

The quantities jj.^ j and Gq \ in eq.s d25l > and d26| > can be computed for any glitch instance 
in the network, i.e. for any allowed set of (possibly null) SG atoms in the interferometers' 
outputs. For each instance, these quantities identify the corresponding marginal moments 
of the detection statistics C,- under Hq. We are obviously interested in computing the 
same moments averaged over all possible glitch realizations in the network detectors, using 
the known prior distributions for the glitch parameters. After some tedious algebra, we 
accordingly get (under the usual discrete-time representation): 

Moo = °> 

> 2 (27) 



(og) =f- 2 (Aj+A})c*N s ., 



4=0 



{4 



l 

4l = -fr 1 (Aj+A e )E( W . ¥ ' T ), 
(o$) 2 - 3 (tr ( J) 2 +f s 7 2 a 2 [a{A) +Aj)E( ¥ - W T )+ 2AjA e E{ W ■ W ' T ) 
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(') 

^0,3 



-f-\A j+ A e )E( V - V ' T ) 



(°$) =(<$) +fr 2 <T 2 [2(A 2 j +A 2 e )E( W - W T )+2A J A ( ,E( W -y> T ) 



(30) 



where y — {y/(?i), i//(? 2 ) , • • • , W{ l N s )}, and y T ) and £■(<//• </ r ) (multiplied by /7 1 ) are 
the expected glitch energy and the expected correlation between glitches occurring in different 
detectors, respectively, both expectations being taken over all possible glitch instances. 

Note that eqs. (|27| > are nothing but the moments of C, under Hq in the absence of glitches, 
i.e., due to the Gaussian noise floor only. It is therefore apparent that glitches have a twofold 
effect, making the expected value of C; non-zero, and increasing its variance. 

Substituting eqs. d27]i-(l30ll into eqs. (l25"l l and (1261 . the (marginalized) first two moments 
of the detection statistic under Hq can be written (z = 1,2,3) 



and 



where 



Mo 



-U 2 f- l (A j+ A e )E( V - V ' T ), 



1AjA t 



(Af 



-4) 



n 



(31) 



,(32) 



(33) 



4. The Impulsive Component - Toward a Rigorous Approach 

In this Section we relax the assumption made in the previous Section that no more 
than a single glitch may occur in the analysis window in each interferometer, and adopt a 
general, fully rigorous approach to model the glitchy component, along the lines laid out by 
D. Middleton in a series of seminal papers ifTTl . 

The impulsive noise component in each interferometer is accordingly modeled as a 
random process consisting of a linear superposition of atoms, vizQ 

8i(t)= £V(«-$ ) : a i* ) )'' e0 '' »'= 1,2,3. (34) 
k=i v 7 

Here is the chosen representation atom, t^} are a set of random glitch firing times, af^ is 
a set of random (independent) shape-parameters (e.g., amplitude, center frequency, duration, 
bandwidth), and K[[T] is also a random variable, denoting the number of glitches occurring 
in the analysis window ©,, whose time-width is denoted as T. The firing-times and shape- 
parameters are determined independently at each glitch occurrence (i.e., for each k). 

The key modeling assumption is that the glitching component may be taken as stationary 
(homogeneous) on time scales sufficiently long compared to the analysis window. On such 
time scales typical glitches will show up with constant probabilities, and occur at a constant 
rate, which can be stimated from actual data. 

The number of events Ki[T] will be accordingly ruled [28| by a Poisson distribution, i.e. 

prob {K, [T] =Ki} = -4— , (35) 

A;! 

J A straightforward generalization is obtained by adding several terms like (34) using different atom families. The 
characteristic function of the resulting process will be the product of the characteristic functions of its terms. 
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where Nj (aka, A,-7\ A,- being the glitch firing-rate) is the average number of glitches occurring 
in interferometer-/ in the (T-seconds wide) analysis window. We shall assume the above 
probability, as well as the distributions of the firing-times and shape-parameters in d34l i, to be 
the same for all instruments, and henceforth drop the index i. 

The characteristic functions of the process d34l i can be computed exactly up to any order 
IfTTl . The first order one can be written 



F g (S,t) = £ prober] = K}F g ^,t\K), 

K=0 



(36) 



where F g , t \K) is the conditional characteristic function, given K glitches in the analysis 
window ©, viz.: 



F g £,t\K)=E{exp 



m—l 



(37) 



The expectation in d37l > is taken with respect to both the firing times, fl , and the shape 
parameters, cf m \ The pertinent distributions being assumed as time-invariant in 0, and 
independent for each glitch occurrence, eq. d37b and (|36"1 > become, respectively 



F g (^,t\K)=E 



and 



F g ^,t) = exp[N^E 



- 1 



(38) 



(39) 



From the characteristic function F g (§ ,t) it is straightforward to compute the moments of 
the process g(t), representing the impulsive (glitch) noise component in each interferometer: 



Mi 



(Q) 



yielding 



and 



Mi 



(i) 



nP=E[g(tf] =N^E z [\i/{t-t ;a)}+NE[\if z (t-t Q ;a)} 



■■E[g(t)]=NEMt-to;ct)}, 

Tt 2 t,2 



(40) 



(41) 



(42) 



where, the expectations are taken with respect to both to and a. The related distributions being 
assumed as time-invariant in 0, the moments (RTt . d42l) are also time-independent. 



4.1. RK Correlators Distributions under Hq - Rigorous Approach 

Using the model exploited in Section H] for the impulsive component of the instrument 
noise, it is possible to compute the first two moments of the distribution of C, under Hq in a 
rigorous way. Formally, these are obtained by making the substitution 

n(f) —►«(»)+*(*), (43) 

for the noise in each detector in computing Zs[C,- |i/o] an d V ar [Ci\Ho], thus obtaining (to second 
order in the noise moments of g) 

Mo (0 = E[Q\H ] = -f- l N s (Aj +A e ) (44) 
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and 



(a«) 2 = Var[Q\Ho} = (A 2 +A 2 )o*f- 2 N s 




(45) 



and then using eq.s (|4TT > and d42l for the first two moments of the impulsive components gj(t ) 
in (|45j to get: 



$ =E[Q\Ho] = -f- l N s {Aj+A t )N 2 E 2 [y{t~t ;a)l 



and 



l)\2 



■ Var[C i \H ) = (a^ 0) y 



[\+,yf{D. s )]N 2 



N 



E[\(f 2 (t-to;a)} 



y(t-to;a) 



(46) 



(47) 



It is now interesting to compare eq.s d46b and d47l i to eq.s ( f3Tb and d32l ). obtained from the 
heuristic reasoning in the previous Section. In order to do so, the sum in (l36b should include 
only the K = 0,1 terms, to match the assumption made there that no more than a single glitch 
may occur in the analysis window. This gives the following approximate expressions for the 
first two moments of the impulsive noise component, 



E[g(t)\K = 0,1] =Ne- N E[y{t~t Q \a)], 
E[g 2 (t)\K = 0,1]= Ne-"E [ W 2 (t - 1 ; a)] , 
yielding, upon substitution in d44l) and d45l ). 



and 



(48) 
(49) 



E[Q\H ,K = 0,1] = -f^NMj +A e )N 2 e- 2N E 2 [ W (t - t ;a)], (50) 



Var[Ci\H ,K = 0,1] = ( 



(0,0) 



1+2 



Ne- N E[y 2 (t-to;a) 



N 2 e- 2N E 2 [\lf(t-t ;a) 



Equations (l50l . < T5TT > reproduce eqs. ( TJIT l. d32l iff 

YlN- l E{\i/-V/ T ) =Ne- 7J E[\i/ 2 {t-t ;a)], 
ifN'^Ei \\r-\\r lT )= ]vV 2F £ 2 [\j/(t -t ;a)]. 



+ 



(51) 



(52) 



Both equalities are trivially proven, noting that for Poissonian distributions II = Ne~ 

In conclusion, the rigorous approach sketched above agrees, in the appropriate limit, with 
the result obtained in the previous Section from a simple heuristic argument. The rigorous 
approach, on the other hand, allows any number of glitches in the analysis window in each 
interferometer, in a natural way. 
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For the special case of SG atoms, assuming a uniform distribution of the glitch firing 
time over the analysis window, we have 



,(1) 



Mr = 0, /if = 

: re: 

against Monte Carlo simulations. 



■ t (l-e- 2 * 2 f$°? 



(53) 



8 IsflT* 

The analytic results for the moments in eqs. d44b . d45t . d53l were checked successfully 



5. RK Correlator Based Detector Performance 

Under the made assumption of Gaussianity of the distributions of the C; under both H\ 
and Hq, it is straightforward to obtain the Receiver Operating Characteristics (ROCs) |19|, 
which completely characterize the RK-correlator based detector. 

In the appropriate surveillance context, the detection thresholds are determined, 
according to the Neyman-Pearson criterion, from the prescribed false alarm probability a 
as 



r 



4 ) erf C -\a) + 4 ) , '=1,2,3, 



(54) 



Note that, in view of eqs. d46*b and d47l i. the thresholds depend on the variance of the Gaussian 
noise floor, the DOA, the number of samples in the analysis window, and, in view of eqs. d4l1 >. 
d42b . the average glitch energy and firing rate. 

The corresponding false dismissal probabilities are 



/i® = 1-P$' = 1 -erfc 



(0 



i= 1,2,3, 



Combining equations d54i i and ( f55l l, we get the explicit expression of the ROCs 



P% = erfc 



H Mi 



(0 



, i = 1,2,3 



(55) 



(56) 



where P^ is the detection probability. 

In the numerical experiments illustrated below we consider the three-detectors network 
including LH, LL and V, using these subscripts accordingly. 

Equations d54l i and (T55T > have been used to obtain the ROCs (/3 vs a curve shown in 
Figures 3 to 8. All figures refer to the + polarized case. Similar results are obtained for the x 
polarized case, and are not reported for brevity. 

It is interesting first to illustrate how the mere noisiness of the pseudo-templates (O spoils 
the performance of the RK detector compared to the perfect matched filter. To do so, we shall 
momentarily ignore the impulsive noise component, by letting N = 0. 

Figure 3 shows the ROCs of the RK correlators (best DOA assumed) for different values 
of the intrinsic SNR 5/, of the GWB. The ROC for the perfect matched filter corresponding 
to SNR = 7, which is conventionally considered as the lowest operational value for this latter, 
is also shown. The value 5/, = 7 corresponds, in the pertinent best DOAs, to 8s =6.22 (LH), 
6.48 (LL), 6.68 (V). 

It is seen that the RK detector's performance is better when using Clh or Cll as a 
detection statistic. In this case, at a given false-alarm rate, one needs roughly to double the 



tt Strictly speaking, these are not ROCs, according to the usual definition, although ROCs can be trivially derived 
from them. 



Modeling the Impulsive Noise Component.. 



12 



SNR to obtain the same false-dismissal level as the perfect matched filter. Using Cy, instead, 
the SNR must be larger by a factor ~ 4. This is due to the different directional response of 
Virgo compared to the two LIGOs, due to its different orientation. This is illustrated in Table 
I, where we collected the values of the quantities 

m si' 1 |r (0| ,0 4" m\ (57) 

representing the signal to noise ratios, normalized to the intrinsic SNR of the incoming GWB, 
of the data (p) and the noisy template (pY), respectively. Table-I shows that in the DOA 
ranges where Virgo exhibits the largest response (largest normalized SNR p), the two LIGOs 
respond poorly, and the pseudo-template obtained from them has a low normalized SNR pr- 
Conversely, in the DOA ranges where either of the LIGOs has the largest response (largest p), 
the pseudo-template constructed from the other LIGO and Virgo has still a decent normalized 
SNR p T . 

Figure 4 shows the ROCs for Clh for 8/, = 10 and three typical durations of the analysis 
window (20, 40 and 100 ms). The ROCs for Cll and Cv are similar, and are not shown for 
brevity. 

Figure 5 is the same as Figure 3, except that here the performances are averaged over the 
whole celestial sphere. In this case, the RK correlator based detector performs worse than the 
perfect matched filter roughly by a factor of 3 in terms of SNR. 

The effect of instrumental glitches is illustrated in Figures 7 and 8. 

In order to draw these figures, we estimated the parameter distributions of the SG-atoms 
to be used in (f34t from (unclustered) triggers collected in 1 week of S5 data, kindly provided 
by S. Chatterji [29 1 . The distributions obtained for fy and a t are sketched in Figure 6. The 
SG-atom amplitude distribution was assumed as uniform in an interval set by the maximum 
SNR in each detector, beyond which the data are vetoed-out, denoted as SNR veU) . 

In Figure 7, the ROCs for Clh, Cll, and Cv are shown for different values of the 
glitch firing rate X, and compared to the no-glitch case. Obviously, as X increases, the 
best achievable false-dismissal vs. false-alarm probability trade-off deteriorates. Here the 
glitch amplitude si assumed as being uniformly distributed up to level corresponding to 
SNR vetu = 100. 

Finally, in Figure 8 the way the chosen SNR veto value affects the performance is 
illustrated. 



6. Conclusions and Directions for Future Work 



We modeled the impulsive noise component following Middleton, using an atomic 
representation for the glitch population. The proposed model allows to describe analytically 
the detector's performance in the presence of glitches. 

Based on the above, we also presented a simple, fully analytic characterization of the 
RK coherent network data analysis algorithm for detecting unmodeled GWBs with known 
DOA and TOO. Under a reasonable central-limit hypothesis for the RK detection statistics 
distributions, we derived and discussed the detector's performance, in terms of its operating 
characteristics. 

Our main results can be summarized as follows. The presence of noise in the pseudo- 
templates spoils the deflection, compared to a perfect matched filter. The related degradation 
factor depends on the direction of arrival, the energy of the signal, and the length of the 
analysis window. The detection threshold, on the other hand, depends on the variance of the 
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Gaussian noise floor, the DOA, the number of samples in the analysis window, the average 
glitch firing rate, and the maximum allowed (veto dependent) glitch energy. Constant False 
Alarm Rate (CFAR) operation is possible, and the RK detector turns out to be reasonably 
robust against instrumental/environmental glitches. 

More or less straightforward developments of this work include i) using a better detection 
statistic, e.g., a linear combination of the RK correlators with (DOA-dependent) coefficients 
chosen so as to maximize the deflection, and ii) allowing for a time-varying glitch firing-rate 
(Cox processes [30|). 

As possible directions for future work we mention i) identifying a better atom dictionary, 
and characterising more accurately the prior distributions of the relevant parameters using 
a systematic matching-pursuit based analysis ll25l of the available glitch databases, and ii) 
exploiting in full Middleton's model to derive more efficient implementations of the detector. 
In this connection, we note that the straightforward extension of the matched correlator to 
unmodeled waveforms provided by the RK algorithm is likely to be not optimal in view of the 
non-Gaussian nature of the instruments noise, whereby some suitable pre-conditioning of the 
data will be most likely required ll3D . 

Finally, we mention the possibility of integrating Middleton's model in a full-fledged 
interferometer noise simulator including glitches. Work along these directions is in progress. 
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Figure la - Sky maps of the factor E, (eq. d23l l) for the RK correlators C, corresponding 
to LH (left), LL (mid) and V (right), for + (top), and x (bottom) linear polarizations. T = 40 
ms; f s = 4096 Hz; 8h = 10. All panels, X axis : 4 . [rad]; Y axis : 9 S [rad]. 



Figure lb - Sky maps of the factor E; (eq. ( 1231 ) for the RK correlators C,- corresponding 
to LH (left), LL (mid) and V (right), for + (top), and x (bottom) linear polarizations. T = 40 
ms; f s = 4096 Hz; 8 h = 100. All panels, X axis : [rad]; Y axis : [rad]. 



Figure 2 - SG atom with go = 1, to = 0.5 s, fo = 100 Hz, a, — 0.02 s. Top: time domain 
waveform; bottom: time-frequency (Wigner-Ville) representation. 



Figure 3 - Performance in terms of ROCs of the RK pseudo-correlator. Optimal DOAs. 
T = 100 ms, f s = 4096 Hz, several 5/, values. Left: Clh', mid: Cll', right: Cy. 



Figure 4 - Performance in terms of ROCs of the RK pseudo-correlator. Optimal DOAs. 
f s = 4096 Hz; 8 h = 10; several T. C LH only. 



Figure 5 - All-sky averaged performance in terms of ROC curves of the RK pseudo- 
correlator. T = 100 ms; f s = 4096 Hz. Left: C LH \ mid: C LL ; right: Cy. 



Figure 6 - Histograms of SG-atom parameter distributions from 1 week of S5 data l29l ). 
Left: center frequency (/o); right: effective duration (c,). 



Figure 7 - Performance in terms of ROCs of the RK pseudo-correlator for different glitch 
rates. Optimal DOAs. T = 100 ms; f s = 4096 Hz; 5 h = 15; SNR vet0 = 100. Left: C LH ; mid: 
Cll; right: Cy. 



Figure 8 - Performance in terms of ROCs of the RK pseudo-correlator for different 
SNR V eto levels. Optimally oriented source. T = 100»m; f s = 4096 Hz; 8 h = 15; N = 0.1. 
C LH only. 



Table I - The quantities in d57l > evaluated at optimal DOAs. 
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